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' We discuss recent results on color superconductivity in QCD at large chemical 

' potential. 
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. The phase structure of matter at non-zero baryon density has recently attracted a 

great deal of interest. The behavior of hadronic matter in this regime is of interest 
C*") , in connection with the structure of compact astrophysical objects and the physics 

of heavy ion collisions in the regime of maximum baryon density. In addition to 
that, we are addressing a fundamental question about QCD, namely, what is the 
' ultimate fate of hadronic matter as we keep compressing it to higher and higher 

densities. And finally, understanding the structure of matter at finite density is 
part of our quest to understand the full richness of the QCD phase diagram as 
we vary the temperature, the density, and the number of quark flavors and their 
masses. 

, At very high density the natural degrees of freedom are quasiparticles and holes 

in the vicinity of the Fermi surface. Since the Fermi momentum is large, asymp- 
totic freedom implies that the interaction between quasiparticles is weak. In QCD, 
because of the presence of unscreened long range gauge forces, this is not quite true. 
Nevertheless, we believe that this fact does not essentially modify the argument. 
We know from the theory of superconductivity that the Fermi surface is unstable 
in the presence of even an arbitrarily weak attractive interaction. At very large 
density, the attraction is provided by one-gluon exchange between quarks in a color 
anti-symmetric 3 state._QCD at high density is therefore expected to behave as a 
j_i ■ color superconductor li BtJ. 

Color superconductivity is described by an order parameter of the general form 

<P={ijFCY D \cT F ^). (1) 

Here, C is the charge conjugation matrix, and Td> Ac, tf are Dirac, color, and flavor 
matrices. Except in the case of only two colors, the order parameter cannot be a 
color singlet. Color superconductivity is therefore characterized by the breakdown 
of color gauge invariance. As usual, this statement,has to be interpreted with care. 
Local gauge invariance cannot really be broken □. Nevertheless, we can fix the 
gauge and study, in weak coupling perturbation theory, the implications of a non- 
zero vacuum expectation value of the type (|l|). The most important gauge invariant 
consequence of color superconductivity is the appearance of a mass gap, through 
the Anderson-Higgs phenomenon. The formation of a mass gap is of course also 
characteristic of a confined phase. Indeed, it is known that in general Higgs and 
confined phases are continuously connected a. 

In addition to that, color superconductivity can lead to the breakdown of global 
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symmetries. We shall see that in some cases there is a gauge invariant order pa- 
rameter for the breaking of the U(l) of baryon number. This corresponds to true 
superfluidity and the appearance of a massless phonon. We shall also find that for 
Nf > 2 color superconductivity leads to chiral symmetry breaking. 

2 The renormalization group 

In order to assess the relative importance of different instabilities of thr^ quark 
Fermi liquid we employ the renormalization group methods described in Uu. The 
idea is quite simple. The low energy excitations of a degenerate Fermi systems 
are quasiparticles and holes in the vicinity of the Fermi surface, described by an 
effective action of the form 

S eff = Jdtd 3 P V f (id t - (e(p) - e F )) $ + S mt . (2) 

Quasiparticle interactions are represented by a tower of 2n fermion operators. We 
can analyze the general structure of the interactions by studying the evolution of the 
corresponding operators as we successively integrate out modes closer and closer to 
the Fermi surface. The main result of this analysis is that, in general, four fermion, 
six fermion, and higher order interactions are suppressed as we approach the Fermi 
surface. This fixed point corresponds to Landau liquid theory. The only exceptions 
are the well known instabilities in the particle-particle and particle-hole channels. 
In the particle-particle channel the dominant operators correspond to two particles 
from opposite corners of the Fermi surface (pf,—Pf) scattering into (j>'p,~ p'f)- 
This kind of scattering leads to a logarithmic growth of the coupling constant as 
we approach the Fermi surface, and the BCS instability. 

In QCD, a complete set of four fermion operators is given by 

O LL = (^l7o^l) 2 , 0° lr = (^70^X^70^.) (3) 

Ol L = (iPl^l) 2 , O l LR = {iptfipLKipRlipR). 

Each of these operators comes in two color structures, for example color symmetric 
and color anti-symmetric 

(r^ b )(^ d )(6 ab 5 cd ±S ad 6 bc ). (4) 

Nothing essentially new emerges from considering superficially different isospin 
structures, or different Dirac matrices. All such structures can be reduced to linear 
combinations of the basic ones (^), or their parity conjugates, by Fierz rearrange- 
ments. In total, we have to consider eight operators. 

We should emphasize that we have restricted ourselves to massless QCD, and 
local operators invariant under the full SU(3)l x SU(3)r x U(1)a chiral symmetry. 
The-extension to operators that break the anomalous U(1)a symmetry is discussed 
in Era. Also, QCD at high density contains long range interactions due to unscreened 
magnetic gluon exchanges. The inclusion of long range forces is discussed in F°H 11 I 
and in section ^ below. 

The operators (^) are renormalized by quark-quark scattering in the vicinity 
of the Fermi surface. This means that both incoming and outgoing quarks have 
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momenta p\,p2 — ±p and ps,p4 ~ ±<f with \q\ ~ pp. We can take the external 
frequency to be zero. A graph with vertices Ti and r 2 then gives 



GiG 2 / (r^^r 



l)k'k 



(r 2 ), y (r 2 v (5) 



with / = g^r^ 2 l°g(A/_R,/A;7y). Here [Aj\r, A^y] is the range of momenta that was 
integrated out. We will denote the density of states on the Fermi surface by N = 
/i 2 /(27r 2 ) and the logarithm of the scale t = log(A/ r / Auv) ■ The renormalization 
group does not mix LL and LR operators, and it also does not mix different color 
structures. The evolution equations are 

d(G^ L + Gf L ) _ N fnLL , nLL , 2 



dt 3 

d(Gy L - 3Gf L ) _ AT(r ,LL ^ n LL\2 

dt 

d(Gg R + 3Gf R ) 



(G^ + Gf L y, (6) 
N{G^ L -3Gf L ) 2 , (7) 



dt 



0, (8) 
(G% R -Gf R f. (9) 



d(Gy R - Gf R ) _ _ 2N <riL H r>L ,. , 
dt 



In this basis the evolution equations are already diagonal. The coupling G\ = 
Gq L + Gf L evolves as 

Gl ^ = l + W 3) Gl (0)t (10) 

with analogous results for the other operators. Note that the evolution starts at 
t = and moves towards the Fermi surface t — > — oo. If the coupling is attractive 
at the matching scale, Gi(0) > 0, it will grow during the evolution, and reach a 
Landau pole at t c — 3/(NG\(0)). The corresponding energy scale is 

A 7fl = Aavexp(-_JL_). (n) 

This is the standard BCS result. The location of the pole is controlled by the 
initial value of the coupling and the coefficient in the evolution equation. If the 
initial coupling is negative, the coupling decreases during the evolution. The second 
operator in (^J) has the largest coefficient and will reach the Landau pole first, unless 
the initial value is very small or negative. In that case, subdominant operators may 
determine the pairing. 

The particular form of the operators that diagonalize the evolution equations 
can be made more transparent using a Fierz transformation. We find 

O dom = 2(^ L C^ L )($LCi> L ), (12) 
Osub,i = \^lC^ r ){4> r C^ l ) + ..., (13) 

O sub ,2 = ^ l C^l)(^lC^ l ), (14) 

O mar = ^ {iPlCj *Pr) (^rC^O^L ) + ■ ■ ■ , (15) 
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where we have ordered the operators according to the size of the coefficient in the 
evolution equations. We find that the dominant operator corresponds to pairing 
in the scalar diquark channel, while the subdominant operators contain vector di- 
quarks. Note that from the evolution equation alone we cannot decide what the 
preferred color channel is. To decide this question, we must invoke the fact that 
"reasonable" interactions, like one gluon exchange, are attractive in the color anti- 
symmetric repulsive in the color symmetric channel. This can be seen from the 
identity 

{X) ab (X) cd = -(S ab S cd + 5 ad S bc ) - -{S ab S cd - 8 ad 5 bc ). (16) 
3 3 

If the color wave function is anti-symmetric, the Pauli exclusion principle fixes the 
isospin wave function to be anti-symmetric as well. The dominant operator does 
not distinguish between scalar and pseudoscalar diquarks. This degeneracy is lifted 
by operators that break the axial baryon number symmetry. These operators are 
associated with instantons. As we will discuss in more detail in section]^, instantons 
favor scalar diquarks over pseudoscalar diquarks. 



3 The effective potential 

The form of the dominant operator indicates the existence of potential instabilities, 
but does not itself indicate how they are resolved. In this section we wish to 
introduce a simple energy functional that captures the essential dynamics of QCD 
at high baryon density. We shall use this functional in order to analyze the ground 
state of QCD with different numbers of colors and flavors. 

The dominant coupling is a color and flavor ant-symmetric interaction of the 
form 

C = G(S ac S bd - 5 ad S bc ) (S^i - S a S jk ) (V>?C 75 ^) (r k C^f) (17) 

where a, 6, . . . are color indices and i,j... are flavor indices. In the following, we 
shall use the notation ICf^f for the color-flavor structure of the interaction. 

In order to determine the structure of the ground state we have to calculate 
the grand canonical potential of the system for different trial states. Since the 
interaction is attractive in s-wave states, it seems clear that the dominant order 
parameter is an s-wave, too. We then only have to study the color-flavor structure 
of the primary condensate. We assume that the condensate takes the form 

«C7b$) = (18) 

4>f b is a Nf x Nf matrix in flavor space and a N c x N c matrix in color space. 

We calculate the effective potential using the bosonization method. For this 
purpose, we introduce collective fields Af b and Af b with the symmetries of the 
order parameter (|l8|). We add to the fermionic action a term G _1 ICf b ^f Af b Aj^ and 
integrate over the dummy variables Af b and Af b . We then shift the integration 
variables to eliminate the interaction term (|i"7j). So far, no approximations have 
been made. We now assume that the collective fields are slowly varying, and that 
A".' can be replaced by a constant. In this case, we can perform the integration 
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over the fermionic fields and determine the grand canonical potential as a function 
of the gap matrix Afj. 

The integration over the fermions is performed using the Nambu-Gorkov for- 
malism. We introduce a two component field \& = (ip,?p T )- The inverse quark 
propagator takes the form 

The grand canonical potential is now given by 

0(A) = iTr[log(5)] + iA-/C-A. (20) 

In order to evaluate the logarithm, we have to diagonalize the mass matrix A4 = 
JC ■ A. Let us denote the corresponding eigenvalues by 8 p (p = 1, . . . , N c Nf). These 
are the physical gaps of the NfN c fermion species. The trlog(S') term has an ultra- 
violet divergence. This divergence can be_regularized in a way that is consistent 
with the renormalization group approach 113. We find 

".••"(A)=i:{-^('» s (f)-0 + 4 s '}- (21) 

Here, £ is a renormalization scale. The grand potential is independent of £, since 
the scale dependence of the first term is canceled by the scale dependence of the 
coupling constant G. The coupling constant satisfies the renormalization group 
equation discussed above. 



4 Superconductivity in QCD with N c colors and Nj flavors 

From the minima of the effective potential we can now determine the groundstate 
properties in a theory with N c colors and Nf flavors. The simplest case is QCD 
with two colors and two flavors. In this case, we have 

A°f = Ae<%, (N c =N f = 2). (22) 

This order parameter is a color and flavor singlet. This would seem to indicate 
that no symmetries are broken, but this conclusion is not correct. QCD with 
two colors has a particle-anti-particle (Pauli-Gurscy) symmetry which enlarges the 
flavor symmetry of the massless theory from SU(2) x SU(2) to SU(A). Spontaneous 
chiral symmetry breaking by a non-zero quark condensate (ipip) corresponds to the 
symmetry breaking pattern SU(4) — > 5p(4). The coset SU(4)/Sp(4) = S 5 , so 
there are 5 Goldstone bosons. These Goldstone bosons are the three pions, the 
scalar diquark, and its anti-particle. Under the Pauli-Gursey symmetry the chiral 
condensate {^ip) is equivalent to the diquark condensate (e ah ei jipfC 751/} j). In the 
diquark condensed phase, the five Goldstone bosons are the three pions, the sigma, 
and the anti-scalar diquark. 

Both a quark mass term and a chemical potential break the SU (4) symmetry. 
If the chemical potential is zero, but the quark mass is not, then the theory has a 
non-zero chiral condensate. If the chemical potential is increased there is a phase 
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transition at a critical chemical potential fi c ~ m^/2. For fi > [x c , the order 
parameter rotates from the chiral condensate .to. the diquark condensate. QCD 
with two colors is discussed in more detail in 

We would now like to discuss real QCD with three colors. The case of only one 
flavor is special, since in this case a color anti-symmetric condensate cannot be a 
spin singlet. We therefore expect condensation to take place in a spin-1 channel. 
This problem is of some relevance for the behavior of real QCD with two light and 
one intermediate mass flavor. For m s larger than some, critical value, QCD has a 
phase with separate pairing in the ud and s sectors LL3'l£I. Here, we will not discuss 
N f = 1 QCD any further. 

In the case of two flavors, the order parameter 

Ay = Ae 3 <%, (N c = 3, N f = 2), (23) 

breaks color SU(3) -> SU(2). The chiral SU(2) L x SU(2) R symmetry remains 
unbroken. At this level, the Fermi surfaces of the up and down quarks of the third 
color remain intact. Subleading interactions can generate a gap for these states. 
The exact nature of this gap is hard to determine, even in the limit of very large 
chemical potential. 

The most-.interesting case is the situation for three colors and three flavors. As 
shown in li3lL3, these results are relevant to real QCD as long as m 2 s < 4/iA. In this 
case the order parameter takes the form 

A# = A(a?J* - a^), (JV C = JV>=3). (24) 

This is the color-flavor locked phase suggested in 0. Both color and flavor sym- 
metry are completely broken, but there are eight combinations of color and flavor 
symmetries that generate unbroken global symmetries. This is obvious from the 
result We can always undo a flavor rotation by an ST/ (3) matrix U by per- 

forming a subsequent color rotation with the matrix U*. Note that since color 
is a vector-like symmetry, the remaining flavor symmetry is vector-like also. The 
symmetry breaking pattern is 

SU(3) L x SU(3) R x 17(1) v -» SU(S) V . (25) 

This is exactly the same symmetry breaking that QCD exhibits at low density. This 
opens the interesting possibility that in QCD with, three flavors, the low and high 
density phases might be continuously connected ta. We will discuss this phase in 
more detail in section [t]. 

The phase structure of QCD-like theories with more than three flavors was 
studied in L3. The main result is that the general structure found in the case of three 
flavors extends to more than three flavors. The order parameter involves a coupling 
between color and flavor degrees of freedom and the chiral SU(Nf)L x SU(N/)r 
symmetry is broken down to a vector-like subgroup. In general, however, this 
subgroup is not the full flavor SU(Nf)v flavor symmetry that is present at zero 
density. In the case of four flavors, for example, one finds SU(4)l x SU(4)e — ► 
0(i) v . 
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5 Instantons and Color Superconductivity 

The first Quantitative estimate of the superconducting gap was given by Bailin 
and Love □. They assumed that the gap is generated by perturbative one-gluon 
exchange. In this case one has to deal with the problem that the quark-quark 
scattering amplitude is singular at small angles. Bailin and Love assumed that 
in a cold quark liquid both the electric and magnetic interactions are screened at 
some scale A. In this case, the gap is on the order of A ~ /xexp(— 6ir 2 /(g 2 L)), 
where L = log(^ 2 /A 2 ) depends logarithmicly on the screening scale. This result 
implies that A is quite small, unless we are willing to consider very large coupling 
constants, a s = g 2 /(47r) > 1. Bailin and Love concluded that, most likely, A//x ~ 
O(10 -3 ) for chemical potentials fx < 1 GeV. We will see in the next section that 
the correct result for the perturbative gap is A ~ /i<? -5 exp(— 3n 2 / (y/2g)) , which is 
parametrically bigger than the result of Bailin and Love. Nevertheless, it remains 
unclear whether perturbation theory is applicable for chemical potentials that are 
of practical interest. 

After the work of Bailin and Love, there was practically no interest in quark su- 
perconductivity until two groups independently pointed out that non-perturbative 
effects, instantons, can generate gaps that are muxiL larger, on the order of 100 MeV 
at densities a few times nuclear matter density oE3. In general it was realized that 
any kind of effective quark interaction which is consistent with the magnitude of 
the chiral condensate and general features of the hadron spectrum at zero density 
will lead to gaps in that range. As we saw above, this is most obvious in the case 
of two colors, where the diquark condensate is related to the chiral condensate by 
a symmetry. 

Instantons are classical solutions of the euclidean Yang-Mills field equations, 
describing tunneling events between topologically different vacuum configurations. 
Instantons are the strongest non-perturbative fluctuations in the QCD vacuum; 
when lattice gauge configurations are made somewhat more smooth, it is instantons 
one mainly finds, see c3 for a recent review. Instantons have a dramatic effect on 
quarks. In the background field of an instanton, the Dirac operator has an exact, 
chiral, zero mode. This zero mode is connected with the chiral anomaly but it also 
plays an important role in connection with spontaneous chiral symmetry breaking. 
There is evidence that the quark condensate can be understood as a collective state 
built from instanton zero modes. These zero modes, just like the chiral anomaly, 
continue to exist at non-zero chemical potential. In this case, the zero modes lead 
to a pairing interaction, and superconductivity. 

This is most easily discussed in terms of effective interactions. In the case of 
two flavors, the effect of zero modes on the propagation of quarks can be encoded 
in a four fermion interaction 

c = 9q„hq { 2J ^ [i^r-^f + (^t 75T -^ )2] 

+ i ^(^Vr-^>) 2 } > ( 26 ) 

where N c is the number of colors, t~ = (t, i) is anJsospin matrix, and T is a form 
factor that reflects the zero mode wave function E3n3. From this interaction we can 
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Figure 1. Quark constituent mass and superconducting gap from a mean field calculation in 
N f =3 QCD. 

directly read off the effects of instantons in different quark-anti-quark channels. In 
the pseudoscalar channel the interaction combines attraction for isospin 1 (pion) 
channel with repulsion for isospin (j]'). Similarly, one finds attraction in the 
scalar isospin-0 (a) channel (responsible for spontaneous chiral symmetry breaking) 
together with repulsion in the scalar isospin- 1 channel (do). 

Using simple Fierz identities, we can construct the effective interaction for color 
antisymmetric [3] and symmetric [6] diquarks. The result is 



where Ti is the anti-symmetric Pauli matrix, and Xa,s are the anti-symmetric and 
symmetric color generators. In the color [3] channel, the interaction is attractive 
for scalar [ip T Cj^i/)) diquarks, and repulsive for pseudoscalar (ip T Cip) diquarks. 

We pup now use this interaction in order to study the phase structure of Nf = 2 
matter E3£aEJ. This is most easily done in the mean field approximation. At zero 
chemical potential chiral condensation is favored over diquark condensation, simply 
because the corresponding coupling constant is larger. At finite chemical potential 
chiral condensation is suppressed because quark loops are Pauli blocked. Quark- 
quark pairing, on the other hand, is enhanced because the density of states on the 
Fermi surface grows. As a result we find a first order transition from a phase with 
chiral symmetry breaking to a superconducting phase. As an example, we show the 
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(27) 



results of a somewhat more sophisticated calculation in Nf = 3 QCD in figure 
[j]. The instanton density was fixed in order to reproduce a constituent quark mass 
AI = 400 MeV at zero chemical potential. In this case, we find the transition to 
the superconducting phase at a quark chemical potential fi c ~ 300 MeV. In QCD 
with three flavors there is a flavor singlet and a flavor octet gap in the high density 
phase. The singlet gap peaks at around 80 MeV. As we will discuss in more detail 
below, chiral symmetry remains broken in the high density phase. There are two 
types of constituent quark masses in the high density phase, but both of them are 
small, M < 10 MeV. We should stress that the mean field description provides a 
very simple minded picture of the phase transition. The critical chemical potential 
is very small, and there is no nuclear matter phase in which quarks are clustered 
into nucleons. 



6 Superconductivity from perturbative one-gluon exchange 

If the chemical potential is very large, asymptotic freedom implies that the coupling 
is weak. In this case, the gap can be calculated in weak coupling perturbation 
theory. The result A ~ fJ-g 5 cxp(— 3TT 2 /(^/2g)) mentioned above was first derived 
by Son using renormalizatkin group arguments. I±_has since then been rederived in 
a number of works Ej'EjO'Ej. Here we will follow E3 and concentrate on the Nf — 2 
order parameter (|25|). The gap matrix has the form 



A?/(g) = (X 2 r b (r 2 ) l3 C l5 ( A 1 (q )l(l + a ■ q) + A a fa,)i(l -3-$)). C>s) 



This decomposition is useful in order to discuss the weak coupling limit and make 
the gauge invariance of the result manifest. In weak coupling, we can replace 
a ■ q — > 1 and only Ai produces a gap on the Fermi surface. We have kept the 
dependence of the gap on the frequency. We will justify the necessity for this below. 
In condensed matter physics, keeping retardation effects in the gap equation is 
referred to as Eliashberg theory. 

We calculate the gap by solving the Dyapn-Schwinger equation for the quark 
self-energy in the Nambu-Gorkov formalism □ 

£(*) = -* 5 2 J -0^r;s( q )r b u D;l( q - k). (29) 

Here, £(fc) = — (S^ 1 (k) — Sq (k)) is the proper self energy, is the quark-gluon 
vertex and D^l(q — k) is the gluon propagator. The off-diagonal components of 
( p9[ ) determine the gap, while the diagonal components fix the wave function renor- 
malization. Since the fermion momenta are large we will ignore wave function 
renormalization. For the same reason, we can also ignore vertex corrections. 

The gluon momentum, however, can become soft and corrections to the gluon 
propagator have to be taken into account. The gluon propagator in a general 
covariant gauge is 

D,u(q) = 4^r + 4^-^> (30) 
q z — G q — .r 
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where P^ L are transverse and longitudinal projectors and G, F are functions of 
energy and momentum. £ is a gauge parameter that should not appear in any 
physical result. 

We can derive two coupled equations for the gap parameters Ai,2(po)- in the 
weak coupling limit, we find that only terms involving Ai have a singularity on 
the Fermi surface. We can therefore drop A2. The resulting equation for Ai is 
independent of the gauge parameter. We find 

.2 r r / I! 



A (Po) = 7^ / dqo / dco S - ——^^£11———— (31) 



■ cost 



12tt 2 J ™J \l-COS0+(G+{p Q -q Q ) 2 )/{2p 2 ) 



1 1 1 



A(g ) 



l-co S 0+(F+(p Q -q o ) 2 )/(2 t i 2 )J ^q$ + A(q ) 2 

The integral over cosO is dominated by small (9, corresponding to almost collinear 
scattering. It is therefore important to take medium modifications of the gluon 
propagator at small momenta into account. For q — > and to leading order in 
perturbation theory we have 

F = 2m 2 , G =f m2 f> ( 32 ) 

with m 2 = N f g 2 n 2 / (4tt 2 ) . In the longitudinal part, m 2 D = 2m 2 is the familiar 

Debye screening mass. In the transverse part, there is no screening of static modes, 

but non-static modes are modes are dynamically screened due to Landau damping. 

In our case, typical frequencies are on the order of the gap, qo ~ A. This means that 

1/2 

the electric part of the interaction is screened at q^ — whereas the magnetic 
interaction is screened at qu — ( 7r /4 • m^A) 1 / 3 . Asymptotically, qu -C qE, and 
magnetic gluon exchange dominates over electric gluon exchange. 
We can now perform the angular integral and find 

am =4*1** iog (r-^) ft {q i ^ 

18tt 2 7 \\po-qo\J y/q 2 + A(q ) 2 

with b = 256ir 4 (2/Nf) 5 / 2 g~ 5 . We can now see why it was important to keep the 
frequency dependence of the gap. Because the collinear divergence is regulated by 
dynamic screening, the gap equation depends on po even if the frequency is small. 
We can also understand why the gap scales as exp(— c/g). The collinear divergence 
leads to a gap equation with a double-log behavior. Qualitatively 



g 2 

1 



18tt 



log(f)] 2 , (34) 



from which we conclude that A ~ exp(— c/g). The approximation (Q) is not 
sufficiently accurate to determine the correct value of the constant c. For this 
purpose, we have to solve the integral equation This can be done by converting 
the integral equation into a differential equationEx In the weak coupling limit, an 
approximate solution is given by 

A(p ) ~ A sin f — 9 -=- log (—)), Po > A , (35) 



3%/27r \po 
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with Ao = &/Ltexp(— 37r 2 /(\/2g))- The final result for the magnitude of the gap on 
the Fermi surface is 



We should emphasize that, strictly speaking, this result contains only an estimate 
of the pre-exponent. This estimate is obtained by collecting the leading logarithms 
from electric and magnetic gluon exchanges. We also note that the result shows 
that the perturbative calculation is self-consistent. Asymptotically, the factor \x 
in the pre-exponent overwhelms the exponential suppression factor exp(— c/g) and 
A ^> Aqcd- This implies that the magnetic screening scale qu ~ (t^/^A) 1 / 3 3> 
Aqqd and the perturbative result for the gluon polarization function is reliable. 

For chemical potentials /i < 1 GeV, the coupling constant is not small and 
the applicability of perturbation theory is in doubt. It is known that at finite 
temperature the convergence properties of the perturbative expansion are extremely 
poor. LLhas been argued that the situation at non-zero chemical potential might be 
better lj. The superconducting gap, however, is dominated by collinear exchanges. 
For [i < 1 GeV we have qm — Aqc d and the gap is mainly determined by very small 
momenta. If we ignore this problem and calculate the gap at chemical potentials on 
the order of 500 MeV, we find A ~ 100 MeV. It is gratifying to see that this result 
is quite close to the low-density estimate based on instantons or other effective 
interactions. 

7 Color-Flavor Locking and Quark-Hadron duality in QCD with 
three flavors 

In this section we wish to study the color-flavor locked phase in QCD with Nf — 
3 flavors in somewhat more detail. In the limit where all quarks are massless 
the underlying theory has the continuous symmetry group SU(3)c X SU(3)l x 
SU(3) rxU (1)b- In the color-flavor locked phase the symmetry is broken following 
the pattern 



The breaking of local color symmetry implies that all the gluons acquire mass, 
according to the Higgs effect. There are no long-range interactions. At low energy, 
there is no direct signature for the color degree of freedom. Color is, in this sense, 
confined. 

As explained in section ||, chiral symmetry is broken because color, being a 
vector-like interaction, locks left handed flavor rotation to right handed flavor rota- 
tions. This is an unusual mechanism for chiral symmetry breaking. The condensates 
(ipLipL) and (ipRipR) of left and right handed quarks are quite separate. Indeed, un- 
less we include instantons, they are not even phase coherent. Nevertheless, because 
both are locked to color they are thereby locked to one another. 

Superconductivity breaks the C/(l) of baryon number. This might seem discon- 
certing at first, but the total baryon number of a finite sample is still conserved. 
The meaning of baryon number violation is that there is a Goldstone mode, the 
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(36) 



SU(3) C x SU(3) L x SU{3) R x 17(1)* 



SU(3) C+ L+R x Z 2 . 



(37) 



phonon, that corresponds to small fluctuations in the local baryon density. The 
macroscopic consequence of this fact is easy transport of baryon number, or super- 
fluidity. We can construct a gauge invariant order parameter for superfluidity by 
convoluting several of the primary, diquark, condensates: 

X = (^evkflfitt), # = e°% fc $W*. (38) 

We note that this order parameter carries the quantum numbers of a di-lambda, or 
the hypothetical H-dibaryon. 

A diquark condensate breaks baryon number but there is a residual Z2 symme- 
try that corresponds to flipping the sign of the quark fields. The primary condensate 
also leaves an axial Z2 symmetry, but in QCD this symmetry is broken by instan- 
tons. At large chemical potential instantons and other semi-classical fluctuations 
are suppressed, and the axial Z2 symmetry is almost exact. One consequence of this 
result is that the ordinary chiral condensate (^ip) is strongly suppressed, whereas 
its square, ((ipip) 2 ) is n °t small, proportional to the primary diquark condensate 
squared V3. 

It is interesting to consider how the photon couples to color-flavor locked matter. 
The original electromagnetic gauge invariance is broken, but there is a combination 
of the original electromagnetic gauge symmetry and a color transformation which 
leaves the condensate invariant. The original photon 7 couples with strength e to 
the flavor matrix diag(2/3, —1/3, —1/3). There is a diagonal gluon G which couples 
with strength g to the color matrix diag(— 2/3, 1/3, 1/3). Then the combination 

- 91 + eG 

7 = , 9 9 (39) 
Ve 2 + 9 2 

leaves the color-flavor locked condensate invariant. In the color superconducting 
medium it represents the physical photon. This photon couples to electrons with 
coupling egj \J e 2 + g 2 , which is the unit of charge in the color-flavor locked phase. 
The coupling to quarks receives contributions from both the color and the flavor 
charges. It is easy to see that, in units of the electron charge, these contributions 
always add up to ±1 or 0. The same is true for gluons, too. 

We are now in a position to examine the low energy excitations in the color-flavor 
locked phase. The spontaneous breaking of global chiral SU(3)l X SU(3)r brings 
with it an octet of pseudoscalar Nambu-Goldstone bosons. In addition to that there 
is a massless phonon, corresponding to the spontaneous breaking of U(1)b- If the 
density is very high, there is also an anomalously light pseudo-Goldstone boson 
associated with the U(1)a symmetry. 

In addition to the Goldstone modes, there are also excitations directly derived 
from the original quark and gluon modes. The quarks form an octet and a singlet 
of spin 1/2 particles under the residual SU(3)v- There is an energy gap for the 
production of quark pairs. Dynamical calculations show that this gap is larger for 
the singlet than for the octet. This means that the states in the octet are true 
long-lived quasiparticles, since there is nothing for them to decay into. We can 
think of these states as baryons, since baryon number is only conserved modulo 
2/3. The gluon fields form an octet of spin 1 bosons under the SU(3)v symmetry. 
They acquire a mass by the Meissner-Higgs phenomenon. 
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The universal features of the color-flavor locked state: confinement, chiral sym- 
metry breaking down to vector SU(3), and superfluidity, are just what one would 
expect for Nf — 3 QCD at low density. We have now seen that the same is true 
for the spectrum of low lying non-Goldstone excitations: There is an octet and sin- 
glet of baryons, and an octet of vector mesons. This means that there is no sharp 
distinction between the low-density, hadronic, phase and the high-density, quark, 
phase. In particular, there need not be a phase transition separating the high and 
low density phases of Nf = 3 QCD. Furthermore, we have realized and old goal 
of QCD: we have obtained a weak coupling, but non-perturbative, description of 
hadronic degrees of freedom in terms of quarks and gluons. 



8 Chiral density waves 

Throughout this contribution we have assumed that the dominant instability at 
large density corresponds to the formation of particle-particle (or hole-hole) pairs. 
This is a very reasonable assumption, because only this kind of pairing occurs for 
arbitrarily weak coupling and uses the whole Fermi surface coherently. Nevertheless, 
in strong coupling, or if superconductivity is suppressed, other forms of pairing may 
take place. Obvious candidates are the formatiori.oLLaxger clusters (nucleons, four- 
quark states, dibaryons) or particle-hole pairing Ej'EjEj 

Particle-hole pairing is characterized by an order parameter of the form 

{$(x)il>{y)) = exp(ip • (x + y))T,(x - y), (40) 

where \p\ = pf is a vector on the.Eermi surface. This state describes a chiral density 
wave. It was first suggested in as the ground state of QCD at large chemical 
potential and large N c . The color factors for particle-particle and particle-hole 
scattering are 

c = — + 1 (op), c = N "- - 1 (ph), (41) 

suggesting that particle-particle pairing, and superconductivity, is suppressed at 
large N c . Particle-hole pairing is not suppressed, but suffers from the fact that a 
direction is singled out and only a small part of the Fermi surface is used. At large 
N c , however, screening due to fermions is weak and the perturbative one-gluon 
exchange interaction is even more strongly dominated by collinear exchanges. This 
favors the particle-hole instability. 

Screening was ignored in the original work of cB, but is taken into account in the 
more recent investigations EJ'ELl. The main conclusion is that, in weak coupling, the 
chiral density wave instability requires very large N c 3> 3. We have verified this fact 
from numerical solutions of th^chiral density wave gap equation. In perturbative 
QCD, the gap equation reads 



E(p " )=C l 2 i/ d9 "{ l0g ( 1 



A3 



0|| - <7||) 3 + \ m V\V\\ - 3||l 



2 



o lo § 1 + 71 Z^—T -r==f ( 42 ) 



(PH -q\\) 2 +m%J J ^J q 2 
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Figure 2. BCS (solid line) and chiral density wave (dashed line) gaps from a pcrturbativc 
calculation in QCD with N c colors and Nf = 2 flavors. 



This gap equation is very similar to the BCS gap equation (|3^), except that the color 
factor is different, and that the angular integration is cut off at a scale Aj_ = ^/2jlq^. 
This fact reflects the one-dimensional nature of the chiral density wave. 

Numerical results for both the BCS and chiral density wave (CDW) gap are 
shown in figure ||. We find that for moderate N c , in particular in the case N c = 
3, the chiral density wave instability only exists for very small /i and very large 
coupling. In the case of N c = 3, the CDW gap is never bigger than the BCS gap, 
even if the coupling is big. 

These results are not really quantitative. Only for very large N c does the chiral 
density wave instability occur in weak coupling where the perturbative analysis is 
appropriate. Nevertheless, our results illustrate the point that in the real world, 
with N c = 3, superconductivity is the dominant instability even at moderate chem- 
ical potential. Using more realistic interactions in the strong coupling regime is not 
likely to change this situation, since other interactions, like instantons, will typi- 
cally lack the strong collinear enhancement of the one gluon exchange interaction. 
In addition to that, the true ground state is not determined by the size of the gap, 
but by the magnitude of the condensation energy. Since the chiral density wave 
state does not use the whole Fermi surface, the condensation energy is only of the 
order /iS 3 , as compared to [i 2 A 2 in the case of superconductivity. 
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9 Conclusion: The many phases of QCD 



We would like to conclude by summarizing some of the things we have learned 
about the phase structure of QCD-like theories at finite temperature and chemical 
potential. We begin with the case of two massless flavors, figure ||a. If we move along 
the chemical potential axis at temperature T = 0, there is a minimum chemical 
potential required in order to introduce baryons into the system. Since nuclear 
matter is self-bound, this point is a first order transition: The density jumps from 
zero to nuclear matter density. Along the temperature axis, the line of first order 
transitions eventually ends in a critical point: This is the endpoint of the nuclear 
liquid-gas phase transition. If we continue to increase the chemical potential, we 
encounter the various phases of nuclear matter at high density. Many possibilities 
have been discussed in the literature, and we have nothing to add to this discussion. 
At even higher chemical potential, we encounter the transition to quark matter 
and the two flavor quark superconductor. Model calculations suggest that this 
transition is first order. This is also consistent with the fact that, asymptotically, the 
superconductor is of the first kind, and expected to exhibit a first order transition 
to the normal phase at the critical temperature for color superconductivity. In the 
case of two massless flavors, universality arguments suggest, and lattice calculations 
support, the idea that the finite temperature zero chemical potential chiral phase 
transition is second order. In this case, the linn.pf. first order /i ^ transitions 
meets the T ^ transition at a tricritical point E3£3. 

This tricritical point is quite remarkable, because it remains a true critical 
point, even if the quark masses are not zero, figure ^p. A non-zero quark mass 
turns the second order T ^ transition into a smooth crossover, but the first order 
fi 7^ transition persists. While it is hard to predict where exactly the tricritical 
point is located in the phase diagram it may well be possible to settle the question 
experimentally. Heavy ion collisions at relativistic energies produce matter under 
the right conditions and experimental signatures of the tricritical point have been 
suggested E3. 

We have already discussed the phase structure of Nf = 3 QCD with massless 
or light degenerate quarks in section |^. We emphasized that at T = the low 
density hadronic, phase and the high density, quark, phase might be continuously 
connected. On the other hand, there has to be a phase transition that separates 
the color-flavor locked phase from the T = /x = hadronic phase. This is because 
of the presence of a gauge invariant U(l) order parameter that distinguishes the 
two. In the case of Nf — 3 massless flavors the finite temperature phase transition 
is known to be first order. We expect the transition from the superconducting to 
the normal phase at T and large /u to be first order, too. This means that there 
is no tricritical point in figure [|c. 

The phase diagram becomes more complicated if we take into account the effects 
of a finite strange quark mass, figure |^d. Consider increasing the strange quark mass 
in the color-flavor locked. phase. This shifts the Fermi momentum of the strange 
quarks with respect to the light quarks. If m^/(Afj,) is larger than A(m s = 0) pairing 
between strange quarks and light quarks can no longer take place, and there is a 
first order transition to a phase (2SC) with separate pairing in the ud and s sectors 
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. Asymptotically the gap grows with /i and we expect the color-flavor locked 
phase to dominate for any value of m s . Increasing the temperature or reducing 
the chemical potential favors the 2SC phase. Whether in the case of the physical 
value of the strange quark mass there is an interlude of the the 2SC phase along 
the fj, 7^ axis, instead of a direct transition between the CFL phase and nuclear 
matter, current calculations are not sufficiently accurate to decide. We know that 
there is at least one phase transition, because nuclear matter and the color-flavor 
locked phase are distinguished by a gauge invariant U(l) s order parameter. This, 
of course, is based on our belief that nuclear matter is stable with respect to strange 
quark matter or hyperonic matter. Current calculations have also not conclusively 
answered the question whether the transition along the T / axis is a smooth 
crossover, as indicated in figure |^d and favored by some lattice calculations, or 
whether the transition is first order, as would be the case if m s is sufficiently 
small. This is clearly an important question in connection with the existence of the 
tricritical point. 

Thinking about QCD at finite baryon density has taught us many interesting 
lessons about the phase diagram, but there remain many question marks in the 
picture presented here, and much work remains to be done. 
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